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In this paper we present a comprehensive analysis of the coherence phenomenon of two coupled 
dissipative oscillators. The action of a classical driving field on one of the oscillators is also analyzed. 
Master equations are derived for both regimes of weakly and strongly interacting oscillators from 
which interesting results arise concerning the coherence properties of the joint and the reduced 
system states. The strong coupling regime is required to achieve a large frequency shift of the 
oscillator normal modes, making it possible to explore the whole profile of the spectral density of 
the reservoirs. We show how the decoherence process may be controlled by shifting the normal mode 
frequencies to regions of small spectral density of the reservoirs. Different spectral densities of the 
' reservoirs are considered and their effects on the decoherence process are analyzed. For oscillators 

\ with different damping rates, we show that the worse-quality system is improved and vice-versa, a 

■ result which could be useful for quantum state protection. State recurrence and swap dynamics are 

' analyzed as well as their roles in delaying the decoherence process. 
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I. INTRODUCTION 



The process of decoherence of quantum states has long been a central issue in the description of quantum measure- 
ments [1-3]. In recent years, experimental advances in the domain of cavity QED and trapped ions have allowed the 

• decoherence of photon [4] and phonon [5] field states to be probed in more depth, providing insights into the borderline 
\ between the classical and quantum descriptions of the physical world. The decoherence time of a superposition of 
■ coherent states in a cavity field was measured [4] and shown to be in full agreement with theoretical predictions [6,7]. 

In trapped ions systems, the observed damping of Rabi oscillations has motivated a number of articles on the main 

• sources of noise leading to decoherence [8-11]. Such experimental achievements in matter- field interactions have also 
^ i' encouraged a deep dialog between theoretical and experimental physics, resulting in a degree of mastery of funda- 
^ ' , , mental quantum phenomena that may herald a new stage in the technology of communication [12] and computation 

[13]. 

^ , The exploration of the borderline between quantum and classical descriptions of nature [2] has impelled the gen- 
^ eration of superposition states of mesoscopic systems, known as "Schrodinger cat states" [4,5]. Such superpositions 
are irreversibly affected by their surroundings, whose effect is to destroy probability interference (coherence), and 
continuously transformed into statistical mixtures. Thus, the environment plays a key role in the establishment of a 
direct correspondence between quantum and classical dynamics. While the decoherence time of a superposition state 
depends on the amplitude of the field, the relaxation does not, since the model adopted for the relaxation process is 
' amplitude damping, achieved by coupling the systems bilinearly to the degrees of freedom of the reservoir. 

Decoherence and its dependence upon the amplitude of the superposition state is the main obstacle to the imple- 
mentation of a logic network based on quantum gates [14,15]. The dream of quantum communication and computation 
comes up against the nightmare of decoherence mechanisms [16], owing not only to the inevitable action of the sur- 
rounding environment but also to the intrinsic fiuctuations in the interaction parameters required for logic operations 
[8-10]. The need for huge superpositions of qubit states in the practical implementation of logical operations imposes 
the requirements that the quantum systems be totally isolated from the environment and that the interaction param- 
eters involved be tightly controlled. For this reason, investigation of the sources of noise in such promising quantum 
systems is a crucial step towards the realization of a quantum logic processor. There is also a major effort being 
made in present-day research, to discover mechanisms to prevent decoherence occurring in actual physical systems, 
by using parity kicks [17], stroboscopic feedback [18], engineered driving fields [19] or an engineered reservoir [20-22]. 
In this light, the main concern of the present work is to analyze the coherence dynamics and decoherence process in a 
network composed of two coupled dissipative oscillators, which may be field modes in dissipative cavities [23], phonon 
modes of trapped ions [24], phonon modes of surface electrons in fiquid helium [25], etc. Master equations are derived 
for both weakly and strongly interacting oscillators, leading to interesting results concerning the coherence properties 
of the joint and the reduced system states. This work constitutes a first step towards a more comprehensive treatment 
of the decoherence process in multipartite quantum systems. 
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On attempting to extend the work on dccohcrcnce to interacting quantum oscillators coupled to distinct reservoirs, 
one faces the problem of deriving a master equation for different regimes of coupling between oscillators. Leaving to 
one side the difficulty of engineering an arbitrary coupling strength between the oscillators, in the present work we 
analyze not only the weak, but also the strong coupling regime, where the coupling strength between the oscillators is 
near the typical oscillators frequencies. In both regimes we assume that the coupling strength between the oscillators 
is considerably larger than the system damping rates. 

In weak coupling, the coupling strength between the oscillators, labelled £ = 1,2 from here on, is considerably 
smaller than the typical frequencies of either oscillator, and the resulting master equation is as if the two oscillators 
were decoupled and a decay channel, described by the Liouville operator Cf p, can simply be inserted into the master 
equation for each oscillator considered. In that case, assuming both oscillators have the same damping constant, the 
decoherence time for each oscillator is unaffected by the interaction with the other one. However, when the oscillators 
have different damping constants (field modes in cavities with different quality factors, for example), we observe that 
the "good-quality" oscillator gets worse, while the "bad-quality" oscillator gets better, a result which can be employed 
for quantum state protection. 

In the strong coupling regime, we observe that a cross-decay channel £i2p appears, besides the usual system-reservoir 
individual decay channels Cep. This cross-decay channel modifies the decoherence process of both the joint and the 
reduced system state, to an extent depending crucially on the spectral density of the reservoirs. In fact, in the strong 
coupling regime, the normal-mode frequencies are substantially shifted from the typical oscillator frequencies, enabling 
us to explore the whole profile of the spectral densities of the reservoirs. We show how the decoherence process may be 
controlled by shifting the normal- mode frequencies to regions of small spectral density of the reservoirs. Apart from 
these spectral densities, the competition between the cross-decay and the usual channels can give rise to a computed 
delay or advance of the decoherence process, for eigenstates of the system normal modes. Note that if a system pointer 
variable docs not commute with the operator responsible for its coupling to other system, it is clear that the internal 
dynamics must interfere in its decoherence time. Thus the derivation of master equations for strongly interacting 
systems is a central task [26,27] in the study of decoherence in quantum networks. 

It is worth mentioning some previous work concerned with coupled systems. In Ref. [23] the authors describe a 
proposal to achieve reversible decoherence of a mesoscopic superposition of field states. This proposal is based on the 
possibility of performing a reversible coupling between two Fabry-Perot cavities. In Ref. [28], a theoretical model of the 
experimental proposal in Ref. [23] is given, but in [28] the inevitable coupling of the resonators to their environment 
is taken into account when the reversibility of coherence loss is analyzed. A system of two coupled cavities is also 
analyzed in Ref. [29], where just one of the cavities is interacting with a reservoir. In Ref. [29], a master equation is 
derived in the case of strongly coupled cavities and it is shown that the relaxation term is not simply the standard one, 
obtained by neglecting the interaction between the cavities. It is the aim of the present paper, in the context of cavity 
QED, to analyze the reversible decoherence process of Rcfs. [23,28], where two dissipative cavities are considered, 
but investigating also the regime of strongly coupled cavities, as done in [29], in which some remarkable coherence 
properties appear. A central result extracted from our discussion is that in a strongly interacting quantum network 
the decoherence time may not decrease as the number of systems considered grows. In fact, it may happen that with 
many coupled sites [30], as with the two coupled systems analyzed here, the decoherence time increases, depending 
on the spectral density of the reservoirs. 

Together with the strong coupling between two oscillators we consider a classical driving field feeding one of the 
oscillators continuously, which is intended to drive each of the coupled oscillators to a stationary coherent state [19,31], 
other than the vacuum state. We select a specific coupling between the oscillators, which may be responsible for the 
dynamics of local transfer of states or state swap between them, as discussed in Ref. [32] in relation to cavity QED 
with weak coupling. 

The outline of this paper is as follows. In Sec. II we develop a master equation for the two coupled lossy oscillators, 
one of which is under the action of a classical driving field, and we analyze the weak and strong field coupling regimes. 
As the spectral density of the reservoirs plays a crucial role in the strong coupling regime, in Sec. Ill we analyze 
particular cases of spectral densities. In Sec. IV we develop a c-number version of the master equation and solve it 
with reservoir temperatures set to zero. The solution of the master equation derived for strongly coupled oscillators is 
also analyzed in Sec. IV for special initial field states. In Sec. V we discuss the recurrence and swapping dynamics of 
the system states. The central result of the paper, namely the coherence properties of the system states, which depend 
on the spectral density of the reservoir and are strongly affected by the regime of coupling, is presented in Sec. VI. In 
Sec. VII we consider the two oscillators to have different damping rates and demonstrate that the coupling between 
them, assumed to be larger than these damping rates, makes the good-quality oscillator worse and the bad-quality 
oscillator better. A careful analysis of the entropy excess in our network is developed in Sec. VIII. Finally, Sec. IX 
concludes the paper. 
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II. THE PROBLEM OF COUPLED DISSIPATIVE OSCILLATORS: DERIVATION OF THE MASTER 

EQUATION 



General results can be extracted from specific examples of qiiantum oscillators, such as field modes in coupled 
cavities, trapped ions, or surface electrons in liquid helium. Let us consider a system of two interacting oscillators 
under the action of a driving field, as pictured in Fig. 1. We start from a positive-defined Hamiltonian so that the 
energy spectrum has a lower bound which is equal to zero in the absence of the driving field [33,34]. The system 

Hamiltonian is then given by 



A 

-ai 



+hF (e^' + aj) [e-^^ + aj) , (1) 



a2 

l^2k 



where a\ and ai are, respectively, the creation and annihilation operators for the oscillator mode of frequency uj(,o, 
whereas b^k and fej^, are the analogous operators for the fcth bath mode of oscillator whose corresponding frequency 
and coupling strength are uj^k and Viki respectively. The coupling strength between the oscillators is A and the 
classical driving field applied to oscillator 2 has intensity F and frequency uj. Assuming that the coupling between 
the oscillators and their reservoirs satisfies the condition {Ve,k) l^ih ^ '^ia-, and shifting the origin of the energy 
scale to hF, we obtain from Eq. (1) the Hamiltonian H = + Hj, given by 

He = tkJealae + h'^uj(,kb\,.btk + h'^Vtk {aih\^ + a\he,]^ , (2a) 

k k 

Hi = hX (aial + a\a2^ ■ (2b) 
Here, ue is related to the natural frequency coeo by 

(A^ F \ 
1 + -. h 6t2 ) , (3) 
4a;ioW2o W20 / 

such that in the limit of weak coupling between the oscillators {X/ojeo <C 1) and weak amplification process {F/uj20 ^ 1) 
we obtain the natural frequencies cj^ = oj/iq. In this limit, we could have started from Hamiltonian H given by Eqs. 
(2a) and (2b) instead of (1). Note that since we are assuming weak couplings between the oscillators and their 
reservoirs, it is unnecessary to write these interactions in a positive-defined form, as done in Eq. (1). However, the 
positive-defined form for the coupling between the oscillators ensures an energy spectrum with a lower bound (equal 
to zero when the driving field is switched off), whatever the value of the coupling A. Under a unitary transformation 
represented by the operator 



U{t) = exp 



iut alat + J2 , (4) 

^ \ k J. 



we obtain the transformed time-independent Hamiltonian 

n = U\t)HU{t) + i^^^U{t) = Y,^i+ W/, (5) 
dt ^ 

where 

He = hjj'ea\ae + h'^uj^^blifitk + /lA {aia\ + a\a2^ , (6a) 

k 

Hi = Vifc (aebl^ + albek) + hF {a\ + 02) fe, (6b) 
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and the shifted frequencies are given by 



■ oji — u), (7a) 

■ LUgk - L^- (7b) 

From here on we consider the specific case where oj[ = ijj'2 = ^ (or u>i = ^2), which links the ampUtude of the 
driving field with the natural frequencies, as F = (wio — W20) [l + A^/ (4wioa;2o)] • Therefore, in the absence of the 
driving field, the condition uj[ = Lu'2 = fl, implies that cjio = uj2a- With the condition ui'i = (j'2 = O, the Hamiltonian 
in Eq. (6a) can be diagonalized through the canonical transformation 

Ai = ^ (ai + 02) , (8a) 
v2 

-42 = (oi - 02) , 

where Ai and A2 satisiy the same commutation relation as ai and 02: [^i, Aj] = and 



A,A\ 



(8b) 



The purpose of 



these new operators is to decouple the direct interaction between oscillators 1 and 2, described by hX ^0102 + a\a2 

Consequently, indirect interactions between oscillators 1 and 2 will be created through their respective reservoirs, as 
described by Hamiltonian H = Hq + H/, where 



Hn 



^tA\Ai - (-1)'^ (4 + ^Y.^'tk^\M 



= -i E (^^^Ife + - (-l)'^2fe {A,h\^ + A 



\b2k 



(9a) 
(9b) 



and the frequencies arc given by Jl^ = 17 — (— 1)^A. Without direct coupling between oscillators 1 and 2, as modelled 
by Hamiltonian H, it becomes simpler to derive the master equation, following the reasoning in Ref. [26]. In the 
interaction picture, to the second order of perturbation, the evolution of the density matrix of the coupled oscillators 
is given by 

(10) 



dt 



^ dt'TTR [V{t), [Vit'),pRiO) ® pi2{t)]] , 



where V(t) ~ exp (iHoi/ft) H/ exp (— iHot/fi). Note that the density matrix in the interaction picture, px^ify^ follows 
from the state vector transformed by both unitary operators: U'^it) defined by Eq. (4) and exp (zHot//l). Defining the 
operator 0\[t) = 5I^/,/i4fe?4fe ^^P /v^, we proceed to obtain the master equation, assuming that the reservoir 

frequencies are very closely spaced, to allow a continuum summation. We have to solve the integrals appearing in Eq. 
(10), related to correlation functions of the form 

dt' (o\{t)Oi{t')) e-^^^t+i^nt' = dt' 



'^^Kvl {iOik) (t^£fc) [Ni [ujtk) + 1) 



47r 



(11) 



1,2, the fimction Ni^ {i^ek) is defined by 



where, from here on, 



and ai {oJtk) is the density of states of reservoir t. Performing the variable transformations r 
u)e,k — — w, we obtain, for n — 1 and 2, respectively. 



(12) 



t — t' and e 



^ dt'(^Ol{t)Oe{t'))e- 



1 



= - e 



/oo 
d£y/(. + a;/)a|(£ + a;+) 

X [Ne (e + w+) + 1] / 

JO 



dre" 



j\t' {0\{t)0i{t') 



1 

— c 

4 



i(p.m.-^-2)t 



/OO 

X [Ni {e + uj^) + 1] f dTe''\ 
Jo 



(13a) 



(13b) 
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where the frequency lo^ has been split into two effective frequencies corresponding to the normal modes of the coupled 
oscillators (note that wi = a;2): 



(14) 



We have suppressed the contribution of the Cauchy principal value since it represents only a small shift in the frequency 
w^. We note that the minimmn value of ujj is F [wio/ (wio + W20)] which follows from A = 2uj2o- Next, we discuss 
both regimes: (z) the weak coupling regime \/uJm ^ 1, and {ii) the strong coupling regime A/w^o ~ 1- For the strong 
coupling regime we will assume A/w^o = 2 in order to minimize w^, which becomes zero when the driving field is 
switched off. 

As usual, we consider that V^, ai, and Ni are functions that vary slowly around the frequency u)f, an assumption 

which docs not apply to the function (^7) ^ ["^-^P 

(Huj^/kT) - 1] (taking the reservoir to be in thermal equi- 
librium at temperature T) in the strong coupling regime when F = 0, since in this case u>^ « 0. However, this regime 
(even with F = 0) can safely be applied to a reservoir at absolute zero, the situation we analyze in the present work. 
We observe that, in practice, Ni {t^^) ~ whenever the shift in the frequency ujj , arising from the contribution of the 
Cauchy principal value, becomes sufficiently greater than kT/h. Note that the last integrals in Eqs. (13a) and (13b) 
contribute significantly only when \et\ ^ 1, so that if we extend the upper limit of the time integration to infinity, the 
expressions for the correlation functions become 



dt' lo\{t)Oi{t')) e-^"".*+'^i*' = iMl [N, (c,+) + 1] e-'C^--"!)*, 



^ dt' {o\{t)Oe{t') 
where the damping rates are defined as 



(15a) 
(15b) 

(16) 



Defining = (w^o) o"! (oJio), in the weak coupling regime, where u)f ss « uj^q, we obtain from Eq. (16) the 
result 



1 r 

while in the strong coupling regime, where a;/ ^ , we have 



(17) 



(^7 ) « ^ V^/ (w7 ) g1 (w7 ) when F = 



a) (u>'^ ) when F ^Q. 



(18a) 
(18b) 
(18c) 



^From the above results for the correlation functions, we observe that the master equation for the strong coupling 
regime includes that for the weak coupling regime. In fact, with 7^(0;^) k. and N^{ujf) k. N^{(jJio) we get 

the master equation for the weak coupling regime from that for the strong coupling regime which, described in the 
Schrodinger picture via the mode operators and a\, reads 



dpi2 
dt 



E 



P12, 



A ^ aja„ 



F al 



a2 



+ 



F 



{{ilSii - Xde2) [7<! (i^e ) - It + (^fe - A(5^i) 



2(f^2 _ A2) 
X [27^ (w) - 7^ (^7) - 7^ (w/)] } 
] 

+ 



-\ hi H) + 7« i^e ) {^c )] ( 



a|,pi2 



+ 
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+1 b({^i)+^i{^t )] ( 



aepi2,a 



+ 



It 



K")] ( 



amPi2,aJ 



ae,pi2a 
+ 



+ 



(19) 



We stress that when the driving field is switched off, the second term in Eq. (19) under the summation on £ disappears. 
Otherwise, noting that — = (w^ ~ '^)(w/ — w), when F ^ and cj = , this term becomes 



4A 



Pi2,ae - a 



(20) 



Next, wc derive the master equations in the weak and strong coupling regime. For simplicity, we define 
for the strong coupling regime Jeiojf) = 7^. Besides, we assume from here on the resonance condition for the driving 
field CO = u)2o- 



A. Weak coupling regime 



In the weak coupling regime, where ojf Ki loc k, lo^q, the resonant condition for the driving field implies 7^ (lo) 



7£ 



^ ~ r^/2, so that the master equation becomes 



dpi2 
dt 



= ^ [Pl2,Ho] + ^£^Pl2, 



(21) 



where 



(22) 



and the Liouville operator £^^12 is given by the usual operator structure 



+ 



{Ne M + 1) ( 



aepi2,al 



al,Pi2ae 



ae, Pi2a\ 



)} 



(23) 



Therefore, in the weak coupling regime and assuming that the driving field is resonant with oscillator 2, the interaction 
between the field modes and the amplification process appears only in the von Neumann term of the master equation, 
and does not affect the dissipative mechanism of the individual cavities. However, when the driving field is out of 
resonance with oscillator 2, a correction term is added to the Liouville operator Lipi2, given by 



($72 - A2) 



{Q.5n - A5^i) Y)i (w) - 7^ {u)io)] 



Pi2,ae - a 



(24) 



For a strong amplification process {F/uj^o ~ 1), this correction becomes FQ./ (fi^ — A^) « 1/2 for mode 1 = 2, being 
ignored for mode i=\, since FA/ (9? - A^) < 1. 



B. Strong coupling regime 

Below we present the master equation for the strong coupling regime considering the general situation where both 
driving field and reservoir temperature are present. However, as discussed above, when switching off the driving field 
we must consider reservoirs at absolute zero for our results to be valid. In this regime, where A/w^ ~ 1, the master 
equation is written 

= bi2,Ho] +^£^/9i2 + A2/012, (25) 
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where 



al,Pi2 



, Pi2a| 



al,\pi2,ai\ ^ 



) 



+ 



2 (f22 _ A2) 
X [27^ (a;2o) - 77 - 7^] } 



{(O^a - X6e2) (7^ - 7/) + (f^fe - A^a) 



and a cross-decay channel is included via the Liouville operator 



(26) 



+ [7+Ar,(a;+)-7,-^^4'^7)] ( 



0'mP\2i a, 



[a<!,Pi2aL]) 



4,Pl2 



, ai 



)} 



(27) 



Note that in the weak coupling regime, where « we obtain Eq.(21) from Eq. (25). Evidently, the Liou- 

ville operator accounting for the cross-decay channel, £12/912, owing to the strong coupling between the oscillators, 
can be of the same order of magnitude as the direct-decay channels ll.\p\2 and £2^12- In this regime, in contrast 
to the weak coupling case, a strong driving field (F/a;20 ~ 1) modifies the Liouville operator Lip\2, such that 
F^l (O^ ~ A^) ,-F'A/ (fi^ ~ A^) « 1, independently of the resonance consition. We observe that for both reservoirs 
at absolute zero, the cross-decay channel is lost for the case where 7^ = 7^, which may occur, as discussed below, 
depending on the spectral density of the reservoirs. In what follows it will become clear that the cross-decay channel 
represented by the Liouville operator (27) leads to interesting results concerning the decoherence process in strongly 
coupled oscillators. 



C. The split of the damping rate 

It is interesting to note that for the coupled dissipative oscillators the damping rate 7^ (w^o) = r^/2 for mode I 
splits into 7^ and 7^. To illustrate this mechanism, we assume a Lorentzian coupling Vi between the oscillators and 
their respective reservoirs, such that the damping function Vi (x), centered on frequency xo, is given by 

'^dx) = 'y\-, ^2— ^> (28) 

(X-Xo) +1^ 

with the parameter T accounting for the spectral sharpness around the mode frequency. From the above expression 
and remembering, from Eq. (17), that the frequency a;^ splits into two shifted frequencies w^, we obtain the double 
Lorentzian function depicted in Fig. 2 

r.(x)-^^^^^ 1,1^ 



2 V(x-^/) (x-o;,-) +TV 

= 74'^/)+74'^7)' (29) 

with maxima on w^. In fact, from master equation (19), we observe that in the weak coupling regime, when 7^ (t^/) + 
7^ (w^) = ^ liy^K)) = (since 7^ (i^^) ~ (^7) = 1^*^/2) , we obtain the expected Liouville form for two independent 
dissipative oscillators. From Eq. (29) it is immediately obvious that in the weak coupling regime, where Lof « LOio, 
the damping function presents only one peak, shown by the dotted line in Fig. 2. In this regime, the damping rate, 
assumed to be the maximum of a sharp-peaked damping function, i.e., (for a small value of T), becomes twice the 
value designated for 7^ (w^). As A = ±(w^ — '-^t) increases, the damping function splits into two Lorentzian functions 
whose peak heights arc half the original value F^, as dictated by the master equation (25) and shown by the solid line 
in Fig. 2. The dashed line shows the situation where the two peaks can be clearly distinguished, on the way to the 
strong coupling regime, A/w^o ~ !> where the peak centered on ui^ shifts to around the value F[wio/ ("^lo +'^20)]) 
which can be made as smaller as we wish by decreasing the amplitude of the driving field. In practice, the effect of the 
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strong coupling between the oscillators is essentially to shift the normal- mode frequency lo^ to regions far way from the 
natural frequency of the oscillator lu£o, where the spectral density of the reservoir may be significantly different from 
that around WiQ. In this connection, the spectral density of the reservoir ai (uJik) plays a crucial role in the dissipative 
dynamics of strongly coupled oscillators, since the magnitude of the damping rate {^^e) depends on ae (ive). For 
this reason, we next analyze reservoirs with different spectral densities in order to illustrate the interesting features 
arising from the strong coupling regime. Evidently, the physical systems under consideration and their dissipative 
mechanisms (i.e., the nature of the reservoirs and their spectral densities), will be decisive for our analysis. 



III. SPECTRAL DENSITIES OF THE RESERVOIRS 



It is possible that specific spectral densities could be achieved through engineered reservoirs, a program which 
has recently attracted considerable attention also in attempts to control the decoherence process of quantum states 
[20-22]. Therefore, the results we present below, depending crucially on the spectral density of the reservoir, might 
provide a motivation for future theoretical proposals on engineered reservoirs. For the following discussion we consider 
the strong coupling regime and set both driving field and reservoir temperatures to zero {F, T = 0) and remember, 
from the definition ^i{ujf) = 7^, that the parameter 7^ depends on the reservoir spectral density around ujf. 



A. Markovian white noise 



We start with the simplest case of Markovian white noise, where the spectral density of the reservoir is invariant 
over translation in frequency space, as depicted in Fig 3(a). In this case, assuming a Lorentzian coupling between 
the oscillators and their respective reservoirs, centered aroimd the effective frequencies oj^, as in Eq. (29), we get 
7^ = 7^/2 = r^/4. In fact, the effective frequency shifts to around zero (when F = 0), and so, the system- 
reservoir coupling 7^ becomes half the value 7/ = r^/2, as can be concluded immediately from Eqs. (18a), (18b) 
and 18c. The system-reservoir couplings aroimd uj^ are represented by shaded regions in Fig. 3(a). Therefore, for a 
Markovian white noise reservoir, strong coupling between the oscillators delays the decoherence time of a joint state 
which is an eigenstate of normal mode u>J , as will be discussed latter. Next, we analyze two cases of non-Markovian 
colored noise, still assuming that V^, a^, and are functions that vary slowly around the frequency tof, as discussed 
above. 



B. A Lorentzian spectral density 



Let us consider a Lorentzian spectral density of the reservoir which goes to zero at both effective frequencies u)^ 
and ujf (Fig. 3(b)). The achievement of a reservoir spectral density for which (T£ (w^) ~ would result in a damping 



(x 



+ T2 



0), neither 



(Fig. 3(b)). The achievement of a reservoir spectral density for which (T£ {'-^^^ 
function ji{uif) arising from both terms on the right hand side of Eq. (29) (i.e., Taj/ 

of which would contribute significantly to the relaxation process. Therefore, for such a Lorentzian coupling between 
the field modes and their respective reservoirs, we get from Eqs. (18a), (18b) and 18c the rates ^ F^, since both 
effective frcqiiencics uj^ and shift to regions where the spectral densities of the reservoir are considerable smaller 
than that around ujeo- In view of the; difficulty of obtaining a Lorentzian spectral density sufficiently small around 
w/, such that ae (w/) » 0, we next turn to a more realistic case, which we call a wide Lorentzian spectral density. 



C. A wide Lorentzian spectral density 



Finally, we consider a wide Lorentzian spectral density, which goes to zero at the effective frequency , reaches a 
flat maximum around ujio, continues on it up to around and only reaches zero again at frequencies far beyond ui'^ , 
as depicted in Fig. 3(c). Differently from the Lorentzian spectral density, in this case only the damping function 7^ 

arising from the second term on the right hand side of Eq. (29) (i.e., Tct|/ (x — tUg)'^ + T'^ « 0) does not contribute 

significantly to the relaxation process; i.e., 7^ <^^e, since the resonances around the effective frequency can be 

disregarded. Only the system-reservoir coupling around the effective frequency oj'^ contributes to the relaxation 
process of the cavity mode (7^" = r^/2). Such a model can be considered as similar to subhomic dissipation [3,35] 
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and can be reasonably applied for a variety of phonon-like spectral densities where the de Debye model applies very 
well in the domain of small frequencies. 



IV. THE FOKKER-PLANCK EQUATION 



Using the standard procedures, we derive a c-number version of the master equation (19) for the Glauber P-function 
[36,7] given by 



dP _ 1 ^ 



d d'^ 

+ Q (ai , Q!2 ) 1- Dim a a * 



+ H.C. 



where the function Ce{ai,a2) and the matrix elements Dim satisfy 

Du = it^i (w/) + ijNe, {uj) , 

Di2 = D2i = \Y. bt^i K+) - liNe (w,-)] , 



(30) 



(31a) 
(31b) 

(31c) 



while the parameters 11^, and are defined by 

1 



n^ = 0(7,- + 7/), 

Ei 

Ei 
Bf 



2 
1 

2 

\ (7/ - It + 2iA) , 



iF5t 



F 



2 (f72 - A2 
{n6n - \5t2) (77 - 7/)] 



[(fife - A5a) (27£('^2o) - 7^ -7/) 



(32) 
(33) 
(34) 

(35) 



Note that the drift coefficient Q is linearly dependent upon both variables ai and 02, while the diffusion coefficient 
Dmn is constant, determining an Orstein-Uhlenbeck process [37]. 

In an extreme case of the weak coupling regime, which we do not consider in the present work, where A < F^, the 
drift coefficient C(, depends mainly on ai and the diffusion coefficient is given by Dmn = ^mnDmm, which is strictly 
positive. In this particular case, the dynamics of the two oscillators decouple, since the two-mode state relaxes, due 
to the damping process, before they have time to interact. Thus 



d 

Ill + Ct{at,) — 



D 



u 



dag, 2 dagdal 



+ H.C. 



(36) 



and the general solution is simply the summation of the independent solutions for each mode. 



A. Solution of the Fokker-Planck equation at absolute zero 



^Prom now on we assume both reservoirs set to absolute zero, such that Ng{ujf) = 0. This assumption is equivalent 
to setting Dmn = 0, and thus the Fokker-Planck equation (30) reduces to a simple drift equation 



dP__sr^ 



d 

Ili + Ct{ai,a2)^ + H.c. 



P. 



(37) 



With the substitution P = P' exp [2 (Hi -|- 112) t], we simpliiy the above equation to the form 



dP' 
~dt 



d 

Ce{ai,a2)^^^ +H.c. 



P' 



(38) 
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and assuming that P'{ai,a2,t) = P'{ai{t),a2{t)), we get 

|p'(ai(0,a2W) = I^^A-+H.c^ P'{a,{t),a2{t)). (39) 

Therefore, from Eqs. (38) and (39) we obtain the system of coupled equations 

^ = fii + a^E+ + a^E^ , (40a) 
-^=B2 + a2E+ + aiE^ , (40b) 



which leads to 



^ - {E+ + E+) ^ + ai {E+E+ - E^E^) + {B,E+ - B2E^) = 0. (41) 

Next, we define the parameters (where we have disregarded terms of order 0^(7^/A)) A = 

(7i+ + 72++7r+72")/4 + in, A = {it - it + li ~ I2) /"^^ ^ = (7i+ + 72+ - Ti" - 72") /4, © = 
(7^ — I2 — li +72") /4, and the time-dependent functions 



W = cosh $t l^cos =•= y sin Xtj + i sinh ysin T y cos Xtj , (42a) 
Z"^ = ^sinh cos =•= y cosh $t sin Xt^ + i ^cosh sin T y sinh M cos Xtj . (42b) 
Defining also the parameters 

B2E^ - BiE+ 

-^1 -^2 ^ -^1 -^2 

^1^2" - -^2^1^ . . 

- EtEt-E^E- ^''^ 

we finally obtain the solution of Eqs. (40a, 40b), written in the compact form 

ai (t) =e^'[{al-Gi)W+ + {a^ - G2) Z+] + Gu (45a) 
a2{t) = e^* [{a°2 - G2) W- + (a? - d) Z"] + G2, (45b) 

where indicates ae{t = 0). Prom the solution of Eqs. (45a) and (45b) it follows that the Fokker-Planck equation 
(37) can be solved to give 

P(ai,a2,i) =e2(ni+n2)tp(^^^^2,0)U,^„,(t), (46) 

where P{ai, a2, 0) is the P-function for the initial state. The evolved P-function is then obtained by simply replacing 
at by the evolved parameters given by Eqs. (45a) and (45b). Finally, from Eq. (46) the evolved density operator 
related to (19) follows from 

Pi2{t) = J (fai(fa2P{ai,a2,t)\ai,a2){ai,a2\. (47) 



B. Initial joint states 

Next, we obtain the density operator pi2(i) supposing that the modes 1 and 2 are prepared in a superposition of 
coherent states of the form 

\^)=^f±{\^}j,^3j)±\ph,Ph)), m 
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where J\f± stands for the normahzation factor and the subscripts /, II are related to different coherent states partici- 
pating in the superposition. From Eqs. (47) and (48), the corresponding evolved density operator is given by 

Pl2{t) = ^ Cmn{t) |?m, Cm) Cn I , (49) 

m,7i—1 .1 1 

where the labels <; and ( represent the states of modes 1 and 2, respectively. The expansion coefficients read 

Cmnit) =Nl {±lf-'- {Pi |^l^)-exp(-2Remt)(K-r + |z-r) 

X (Pi |/3^)i-«^P(-2R<^{^>*)(l^^r+l^+r) e'(^i^'n+^'r^i^+^'^i+^iii, (50) 

while the states ^ and C, are given by 

ft = e [{pl - Gi) W- - {Pf - G2) Z+] + Gi, (51a) 
Q = e [{Pl - G2) W+ - {Pi - Gi) Z-] + G2 (51b) 

We define the time-dependent functions 

e = Ini{e-^** [{W^r - Pi*) - {Z-y [P]: - Pi*) 

X [e [GiZ- - G2W+) + G2] } , (52a) 
2^ = e-^^«iA}t |^+ ^z-)* [Pi* {Pi - Pi) - {Pi: - Pi*) Pi] 

- {w+y z- [pi {pi: - pi*) - {pi - Pi) pi*] } , (52b) 
e = im {e -^'^ [{w-y {Pi: Pi*) - {z^y {pi - pi*)] 

X [c {G2Z+ - GiW-) + Gi] } , (52c) 
= e-^^'^W^ {W- {Z+y [{Pi - Pi) Pi - Pi {Pl - Pi)] 

- {W-y Z+ [{Pl - Pi*) Pl - Pl {Pl - Pl)] } . (52d) 
Finally, the reduced density operator can be obtained easily from Eq. (49), being 

p^{t)=Mi y: {±iy-'-{pi\piy--^'-'''''^''^''\'^~\' 

m,n—1 ,1 1 

P2{t)=Afi E {±1)'-'- {pi\ pir-^'-''^'''''''^'-^' 

m,n=I ,11 

X (Pl |/3^)--^--^^^>*)|-"l%^^^^l+^J.^i |C„) (C„| . (54) 



V. STATE RECURRENCE AND SWAP DYNAMICS 



Let us suppose that the joint system 1 + 2 is prepared in the superposition state 

|V'i2)=A4(|a)±|-a))i<^h)2, (55) 

which represents a particular case of Eq. (48), where Pj = —Pjj = a and Pj = Pfj = r}. 

In this section we analyze, in the weak and strong coupling regimes, the effect of dissipation on two phenomena: 
i) the joint-system state recurrence and ii) the state swap between the systems. Considering the dynamics of the 
coupled systems, governed by the Fokker-Planck equation (37), we calculate i) the probability that in each system the 
initial state recurs - the recurrence time for both systems being the same - and ii) the probability of a swapping of 
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states between the systems, i.e., the probabihty of oscillator 1 being in the initial state of oscillator 2, and vice-versa. 
Recurrence is guaranteed whenever a dynamical system is closed and its evolution is limited [38], while state swap 
is only possible for specific interactions. Here the two oscillators are coupled by a bilinear Hamiltonian which, in 
principle, allows state swap to occur between the systems [32]. 

Considering the initial superposition (55), the time-evolved joint state following from the density operator (49) 
reads 

Pl2(i) = ^ Cmn{t)\'i'm,C,m) {^n,C,n\, (56) 
m,n=I ,11 

where 

Cmn{t) = Ml (±1)1-*— exp {-2 |a|' [l - (| |' + 1^" |') exp (-2Re {A} t)\ } 

X exp [i + e^^^ + im {e^± + eifl)] } , (57) 

and the probability of recurrence is given by 
PR{t)=Tri2 [Pi2(t)pi2(0)] 

= Ml Yl C^n{t){'^n\{\a)±\-a)),{{a\±{-a\)\^m){Uv)2iv\Cm)- (58) 

m,n=I ,11 

Analogously to the recurrence probability, the state-swap probability is calculated by swapping the labels of each 
initial state of the two fields, as defined in [32]: 

Ps{t) = Tr[p,2it)pi2mii2)^2il)] 

= Nl C„„(t) (r? Iw) (Cm|(|a) ± I - a))i ((a| ± (-a|) ICm) . (59) 

m.n—1 ,11 

In Fig. 4 we plot the state-swap probability Ps{t) (dotted line) and the recurrence probability PR{t) (solid line) 
against the scaled time \t, taking a = r] = 1 as real parameters. In Figs. 3(a, b and c) we set the driving field 

strength F to zero (so that cjio = C1J20 and the minimum of uj~^ , occurring for A = 2cj20j s-lso becomes zero) and 
assume the absence of dissipation, so as to take these figures as references. In Fig. 4(a) we consider the weak coupling 
regime, assuming A/wio = 2 x 10^^ (a somewhat exaggerated ratio used to visuaHze better the strong oscillations of 
probabilities Pij(t) and Ps{t)). We observe (dotted line), as can be deduced from Eq. (59), that the modes will swap 
their states whenever 

Ats = (2n+l)|, n = 0,1,2,..., (60) 

which means that in the swap time ^5, the state of mode 1 becomes while the state of mode 2 becomes the 

superposition (|a) ± | — a))2- From the solid line and Eq. (58), we observe that the joint-system state recurs whenever 

XtR, = mr, n = 0, 1,2, (61) 

i.e., in the recurrence time tR the mode 1(2) returns to its initial state (|a) ± | — a))-^ (1/3)2)- It is evident that the 
shape of Fig. 4(a) results mainly from the small value of A (compared to fl) which defines an envelope function. 

In Fig. 4(b) an intermediate coupling is assumed, such that X/uio = 1. As can be observed, the recurrence 
dynamics remains the same (solid line), while the swap dynamics begins to be affected by the coupling strength. 
In Fig. 4(c) wo consider the strong couphng regime where A = 2uiiq (wiq = u}2o), such that 0/A = 1 (note that 
fi/A = ojiq/X + A/4a;2o)- As in Fig. 4(a), in Fig. 4(c) the recurrence process remains unchanged, still obeying 
Eq. (61), so that the recurrence time becomes smaller due to the strong coupling parameter A. However, the swap 
dynamics is almost completely lost, the remaining oscillations arising from the nonorthogonality between the states 
(|a) ± I - a))i and |/3)2. 

To understand the behavior of PR{t) and Ps{t) in Figs. 4(b and c), we note that in the absence of dis- 
sipation the recurrence and state-swap probabilities are computed as PR{t) = |('^i2(0) \tlJi2{t))f and Ps{t) = 

2 



(V'i2(0) V'i2(i)l 1(2)^2(1)/ ' respectively. For PR{t), the probability amplitude 

(^-12(0) |V'i2(t)) = Ml ((a kz) (r/ IC/) + exp (2zaG2 sin (At)) {-a (77 |C//) 

± {-a |c/) (r? IC/) ± exp {2iaG2 sin (At)) (a |c//) (r? |C//)) , (62) 
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tends to unity under the conditions CiXii r], <;i ^ a (—a), and c// — > —a (a), which are satisfied when 
cos rii_R cos Aiii = 1. This relation imphes that 

= T = IT' 

where n,m, are integers, both being even or odd. In Fig. 4(b), the intermediate coupUng X/wio = 1 follows from 
n = (4/5)m, such that the joint-system state recurs whenever 

Aifl = 8n7r,n = 0,1,2,..., (64) 

a result which explains the shift in the scaled time Xt observed in Fig. 4(b). 
Concerning the state-swap probability Ps{t), we obtain the expression 

(V'i2(0)|i(2)^2(i) |V'i2(t)) = ATI {{a IC/) {v k/) +exp(2iG2asin(At)) {-a |C//) {v k//) 

± {-a IC/) iv k/> ± exp (2iG2Qsin (Xt)) {a |C//) (r? k//» , (65) 

which tends to unity when rj, Q ^ a {—a), and Qj —a (a). These conditions arc satisfied only when 

sinQts sinAts = — 1 and Gi,G2 « 0. For the special case of -F = 0, as in Fig. 4(b and c), the condition Gi,G2 « is 
automatically satisfied (as can be deduced from Eqs. (43) and (44)) and the modes will swap their states whenever 

(2n-l)7r (2m + l)7r 

= -^T- = 20 ' 

with the additional condition that sin [{n — 1/2) tt] = — sin [{m + 1/2) tt]. Therefore, the maxima of the swap probabil- 
ity are eliminated, as observed in Fig. 4(b) (for the intermediate coupling A/wio = 1) due to the fact that the relation 
2n— 1 = 4 (2m -|- 1) /5 (and so sinOts sin Ats = —1) cannot be satisfied for any pair (n, m). However, owing to the os- 
cillations of Ps{t) within the envelop fimction defined by A(= 4f2/5), the neighborhood of these maxima still survives. 
Differently, in Fig. 4(c) (for the strong coupling regime fi/A = 1), the relation sin fits sin Ats = sin^ Atg ^ —1 implies 
that the single maximum of the swap probability is eliminated (since 0/A = 1, there are no oscillations of Ps{t) 
but the envelope function). A heuristic explanation of the behaviors of Figs. 4(b and c) will be provided below. In 
fact, as discussed below, while the superposition state Af± (|q;) ± | — q;))j^ swaps to oscillator 2, the coherent state [77)2 
does not swap to oscillator 1, even though both return to their respective systems. Evidently, the phase mismatching 
between the coupling parameter A and the field-shifted frequencies ojg = il, (when F = 0), represented by the relation 
sin fits sin Xts = sin^ Xts ^ — 1, is the core of this result. 

In Fig. 4(d), also setting F = 0, we include dissipation, taking both systems with the same damping rate 7^ = r/2, 
where F/wio = 2 x 10~^, and X/iVio = 2 x 10~^ (as in the weak coupling regime) and observing, as expected, the 
relaxation of both probabilities Pnit) (solid line) and Ps{t) (dotted line). To compare the relaxation of probabilities 
Pnit) and Ps{t) in both regimes, in Fig. 4(e) we assume the parameters F = 0, 7^ = 27^7 = r/2 (assuming Markovian 
white noise), F/wio = 2 x 10""^, and A/wio = 2 (as in the strong coupling regime). We observe, comparing Figs. 4(d) 
and 4(e), that in the strong coupling regime the fields recur more frequently, within the relaxation time, than in the 
weak coupling regime. Owing to the strong coupling, in Fig. 4(e) the swap dynamics is almost completely lost, as in 
Fig. 4(c). In Fig. 4(e), a dashed-dotted Hne has been drawn at unity in order to display the slow decay of Pnit). 

Finally, in Fig. 4(f) we again disregard dissipation, but turn on the driving field, taking F/W20 = 1 (i-e-, ^10 = 2a;2o), 
and X/uio = 2 x 10~^ (as in the weak coupling regime), and observe a reduction of the swap probability compared to 
the recurrence probability. (Note that the driving field occurs in the von Neumann term of Eq. (19) even though it 
does not influence the Liouville operator.) This behavior can be explained by the result 

G.-G2 = ^, (67) 

which becomes 1/2 for the parameters considered above, preventing the swap probability from being unity, as required 
by the conditions sin fits sin Aig = —1 and Gi — G2 ~ 0. We have assumed in this figure the ratio lv/loio = 10~^, 
instead of the resonance condition cu/uio = 1/2 adopted above, in order to make clear the ampliflcation effects. ^From 

Eq. (65) it follows immediately that the probability Ps{t) is reduced by the factor exp|— 2[F/(0 — A)]^|, which 

explains the maximum value around 1/2 for the probability Ps{t)- It is worth noting that the expression (67) does 
not diverge for F ^ 0, since in the case where il — X (and switching off the dissipation), Eq. (41) becomes 
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whose solution differs from that leading to Eq. (67). 

It can be argued that the state-swap and recurrence dynamics are consequences of energy transfer between modes, 
which in the presence of dissipation is severely reduced, as the two modes tend to reach energy equilibrium with the 
reservoirs and the driving field. Actually, state-swap and recurrence are properties of information transfer rather than 
energy transfer between systems, even though these quantities are generally correlated. However, it can be shown 
[32] that for the coupling between the modes selected above, even if each mode is kept to constant energy, their 
states can be swapped, showing the independence of the two processes. In the presence of the relaxation process, the 
information transfer between modes is reduced, becouse of the absorption of information by the reservoirs. When the 
systems arc strongly interacting, however, the field states recur many times on the scale of A before the relaxation 
takes place. While not preventing the modes from relaxing to equilibrium with the reservoirs, the recurrence has 
striking consequences for short-time-scale dynamics, such as the dynamics of decoherence discussed below. 



VI. COHERENCE PROPERTIES 



So far we have analyzed the dynamics of strongly-interacting modes 1 and 2 to discuss the recurrence and state swap 
processes. Now we analyze, also in the strong coupling regime, with the driving field switched off {F = 0) and the 
reservoirs at absolute zero (T = 0), the decoherence dynamics of the joint state described by the density operator pi2, 
and of the state of system 1 (2), described by pi — Tt2 pi2 (p2 — Tri P12). In this section we also consider the case of 
identical dissipative systems, 7^ = 7^. The case of different decay rates will be analyzed subsequently. We considc;r 
three different initial joint states: first, the disentangled state, given by Eq. (55), and then the entanglements which 
are eigenstates associated with the normal modes , derived from Eq. (48): with (3} = Pjj = a and /?| = f3}j = —a 
we obtain the eigenstate 

\v>i2) = -Mt {\a, -a) 12 ± \-a, a)^^) , (69) 
associated with the normal mode uj^ , and with /?) = /3| = a and /Sjj = 13} j = —a, we obtain 

^^2) = -^4 (|a, a)i2 ± -a)^^) , (70) 
which is the eigenstate associated with ui^. 



A. Decoherence time of the state IV'12) = A/± (|a) ± | — a))^ ® 177)2 



The coherence of the joint state (55) is given, essentially, by the term 

exp{-2|a|2 [1- (|l^-|%|Z-|^)exp(-2Re{A}t)]}, 



(71) 



coming from the off-diagonal coefficients of the density operator in Eq. (56). In the strong coupling regime, the 
exponential decay in Eq. (71), computed from Eqs. (42a), (42b), reduces to the form 



exp 



2 - e-(^i++^2+)t _ e-(7r+72-)« 



Assuming identical dissipative systems, 7^ =7''', Eq- (72) simplifies to 



exp 



-|a|2(2-e-27+*-e-27-t) 



resulting in a decoherence time for the joint system given by 

TD= [2|a|2(7++7- 



(72) 



(73) 



(74) 



This decoherence time has to be compared with that for an isolated mode in the superposition state J\f± (|a) ± | — a)) 
(also obtained from the weak coupling regime: 7^ = r/2), following from the exponential decay 



exp[-2H2(l-e-r*)] 



(75) 
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and given by the well-known expression « TR/2\a\'^ = Tu, where t/j = stands for the relaxation time of the 
system. Analyzing the decoherence process in the light of the spectral densities considered in Section III, we observe 
that for Markovian white noise (M), where 7" = 7+/2 = r/4, the decoherence time in Eq. (74) rises to 

« 4T^/3. (76) 

For the Lorentzian spectral density (L), where 7^ = ^ F, we obtain the result 

T^«2Tj5/(£+ + £-), (77) 

which becomes large as the reservoir spectral density decreases around the effective frequency uy^ = oj^. We stress 
that for the case where = e~ (i.e., 7+ =7") the cross-decay channel is null. 

Finally, for a wide Lorentzian spectral density (WL), where 7" = e~T/2 <C F and 7+ = F/2, we obtain from Eq. 
(74), the value 

r^^ « 2Td. (78) 

Therefore, for a wide Lorentzian spectral density we obtain a decoherence time for strongly coupled systems which 
is twice as long as in the weak-coupling regime. The mechanism behind these improved decoherence times is that 
in the strong coupling regime (where the natural frequency ujiq of the systems splits into two effective frequencies 
= ± A) the spectral density of the reservoir plays a decisive role in the damping rate, which also splits into two 
Lorentzian functions. When the oscillator effective frequency shifts to regions where the spectral density of the 
reservoir is significantly smaller than that around u^q, the damping rate becomes smaller than its value in the weak 
coupling regime. Despite the spectral densities, the competition between the cross-decay and the usual channels is 
the reason for the computed delay of the decoherence process. 

Still regarding state (55) and the strong coupling regime, focusing on the reduced systems, we obtain from Eqs. 
(53) and (54) the density operators 

p,{t)=Ml Y: (±l)^-''""(/3i|/3^>"^^'*^^'e<^^+«^^^MU(.n|, (79) 

m,n— /,// 

where 



jr(Ai) = exp [-(7+ + 72+)i] + exp [-(7f + 73" )i] + 2 cos(2Ai) exp [-(7+ + j+ + 7f + ^2)^/'^] ■ 
Observe that for identical systems, 7^ =7^, the decoherence of the reduced state in oscillatorl is given by 

1 _ i (e-^f^t ^ g-2^-t + 2cos(2At)e-(T++^")*) 



exp 2|a 
while for the state in oscillator 2 it is 
exp-J -2|q;|2 



|-2|a|2 l-i(e-2'/+*+e-2'^"* + 2cos(2At + 7r)e-('^++'^")*) | 



(81) 



(82) 



(83) 



It is easily shown that the expressions in Eqs. (82) and (83), associated with the decoherence times of the states of 
systems 1 and 2, respectively, oscillate below the curve for the coherence decay of the joint system given by expression 
(73). In Fig. 5, assuming Markovian white noise and setting the fictitious ratio A/F = 5 to make the oscillations 
clear, the dashed and dotted lines refer to the decoherence dynamics of systems 1 and 2, dictated by Eqs. (82) and 
(83); respectively. The solid line represents the coherence decay derived from Eq. (73) and the dashcd-dotted line 
indicates the coherence decay for an isolated mode, computed from Eq. (75). Therefore, for the special case of the 
initially disentangled state given by Eq. (55), the decoherence times of systems 1 and 2 coincide with that of the joint 
system. 
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B. Decoherence time of the state |<Pi2) = A/± {\c't~C')-l2 ^ 



Next, we analyze the decoherence dynamics of the entangled state (69) for the joint system. The decoherence 
process of this joint state, in the strong coupling regime, is given by the exponential decay 



exp 



-4|a|^ (l-e-(Ti 



(84) 



which takes into account only the decay rate of the system-reservoir coupling around the effective frequency to^ . At 
this point it is interesting to identify the contribution of the cross-decay channel to the decoherence process, rewriting 
Eq. (84) as 



exp {-4 \af (1 - exp [- (it + li) t/2] exp [E, (7/ " li) t/2] ) } 



(85) 



In this expression, the term exp [Yli {it ^ T7) ^Z^] comes from the cross-decay channel and it is evident that for 
(7/ — 7^) > 0, its contribution makes the exponential decay slower. Assuming identical dissipative systems, 
7^ = 7~ , the expression (84) simplifies to 



exp 



-4|a|'(l-e-2'^"*) 



(86) 



leading to the decoherence time , which is Td/2 in the weak coupling regime, as expected. However, 

in the strong coupling regime, this decoherence time for an entangled state is equal to To for Markovian white noise 
and is significantly improved for both Lorentzian spectral densities, where we obtain 7d/2£~. 



C. Decoherence time of the state ji^ig) = A/± (|a,Q)j^2 ^ I~'*i~'*)i2) 
The decoherence process of the joint state \^i2)i '^'^ the strong coupling regime, is given by the exponential decay 



exp 



-4|ar 1 



(87) 



depending only on the decay rate of the system-reservoir coupling around the efl[ective frequency w/. Rewriting Eq. 

(87) as 



exp {-4 (1 - exp [- {it + U) V2] exp [- {it ' U) */2] ) } 



(88) 



we identify the term exp [— {it ~ll)^/'A coming from the cross-decay channel. It is evident that for 
E^ (7^ — 7^) > 0, its contribution speeds up the exponential decay, in contrast to the above situation where the 
eigenstate of the normal mode loJ is considered. For 7^ = 7~, Eq. (87) simplifies to 



exp 



-4|a|^(l-e-2^+*) 



(89) 



and the decoherence time 



|2 + 

a\ 7^ 



becomes Tb /2 for the weak coupling regime and also for both spectral densi- 
ties, Markovian white noise and the wide Lorentzian spectral noise, in the strong coupling regime. For the Lorentzian 
spectral density, we obtain 7d/2£+. In the light of the above results we conclude, apart from the dependence of the 
decoherence time upon the reservoir spectral densities, that the eigenstate associated with the normal mode is 
less susceptible to decoherence than that associated with tot ■ 



VII. IMPROVING THE QUALITY FACTOR OF A DISSIPATIVE SYSTEM 

Let us now turn to a central result, which is the effect of the coupling A on the decoherence process for states of 
oscillators with different damping rates Fi and F2. We assume that the state of the whole system is given by Eq. (55) 
and that oscillator 2 has a better quality factor than oscillator 1 (i.e., Fi = F F2), where we have the superposition 
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Af± {\a) ± I — a))-^ which we want to protect. The decoherence of state 1, following from Eq. (79), is then governed 
by the exponential decay 



exp|-2|ap l-i(e-^"* + e-^+*+2cos(2At)e-(^++^"W2) | 
while for the field state 2, following from (80), it is governed by 

exp|-2|a|2 1-1 (e-T"* + e-T^*+2cos(2At + 7r)e-(T^+T")*/' 



(90) 



(91) 



These expressions are to be compared with Eqs. (82) and (83), respectively, where we have assumed the same quality 
factors for both oscillators. The decoherence times obtained from Eqs. (90) and (91) have an upper limit given by 
the exponential decay for the joint system 



exp 



lap 



(92) 



which is derived from Eq. (72). This exponential decay leads to a decoherence time given by 

TB« [|a|M7++7-)]"\ (93) 

which is twice the result found in Eq. (74), for the case of identical dissipative systems {-yf =7^), whatever the 
spectral density of the reservoir. Therefore, even in the weak coupling regime, the decoherence time of a system is 
improved when it is coupled to another system with a better quality factor. We stress that this conclusion holds only 
when A 7^, a situation we have assumed even for the weak coupling regime. Evidently, when A ^ 7^, so that the 
recurrence-swap dynamics does not take place effectively before the relaxation time, the quality factor of a system 
cannot be improved by coupling it to another system of better quality. 

In Fig. 6, assuming the weak coupling regime, 7^ = Te/2, with T1/T2 = 10^ and A/Fi = 5 (a ratio chosen to show 
the dissipative dynamics clearly), the dashed and dotted lines refer to the decoherence dynamics of systems 1 and 
2, derived from Eqs. (90) and (91), respectively. The solid line represents the coherence decay for the joint system 
derived from Eq. (92) and the dashed-dottcd line indicates the coherence decay for the isolated mode 1, where the 
superpositon state J\f± {\a) ± | — a))^ is prepared, computed from Eq. (75). However, while the bad-quality system 
1 gets better when coupled to a good-quality system 2, the latter, in turn, gets worse. In fact, the thick solid line 
in Fig. 6, representing the decoherence process of the superposition Af± {\a) ± | — a))2 prepared in an isolated good- 
quality system 2, displays a slower decay rate than the solid line when system 2 is coupled to a bad-quality system 
1. In conclusion, in a network made of two oscillators with different quality factors, the good-quality oscillator gets 
worse, while the bad-quality oscillator gets better. We hope that this process may be even improved by extending the 
coupling of a bad-quality oscillator to a higher number of good-quality oscillators. 

Finally, we note that the above result can be theoretically employed to protect quantum superposition states 
generated, via atom-field interaction, in open cavities. By coupling an open bad-quality cavity, where we have 
prepared a quantum superposition, to a closed good-quality cavity, we can protect the superposition through the 
coupling between the cavities and, further, in an appropriate recurrence time, we can recover that superposition state 
in the open cavity. The process works as if the original superposition coherence is protected from the dissipative 
mechanism by the recurrence-swap dynamics, much like a system which is put in contact with a high- temperature 
reservoir, but is intermittently brought into contact with a low-temperature reservoir. If the interval between contacts 
is fast enough, the system will take longer to relax to the hot reservoir temperature than it would in the absence of 
the cold reservoir. Of course, this analogy has to be considered cautiously, since temperature and coherence are very 
different features. 

In the next section, we discuss how the original superposition state can be recovered in the bad-quality system after 
its coupling with the good-quality system. 

VIII. ENTROPY EXCESS 

To clarify the role the state-swap processes play in the coherence dynamics, we consider the superposition state in 
Eq. (55), IV'12) = M± {\ct) ± I — ct))i (S) |r?)2, and plot, in Fig. 7(a and b), the probabilities 

Vi{t) = {^j;m\pi{t)\^l:m) , (94) 
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of finding the superposition state = M± (jo;) ± | — 0;))-^ and the coherent state \il>2) = \ri)2 in osciUator t. In these 
figures we assume the strong coupling regime (A/wio = 2), identical dissipative systems (7^ = 7^), and Markovian 
white noise (7+ = 27" = r/2). The ratio A/F = 20 is set to a fictitious scale to make clear the state-swap and state- 
recurrence dynamics, as the strong oscillations obtained with a realistic A would be blurred by the slow (dissipative) 
dynamics. 

^From the solid line in Fig. 7(a), we observe that the superposition recurs to oscillator 1, despite the dissipative 
process. However, the coherent state 177)2 does not swap to oscillator 1, as indicated by the dotted line. In fact, the 
swapping of the coherent state |r?)2 to oscillator 1 would be indicated by the occurrence of maxima of the dotted line 
between those of the solid lino. In Fig. (b) we observe that the coherent state 177)2 recurs to oscillator 2, as indicated 
by the solid line. Moreover, the state |^i) also swaps to oscillator 2 (dotted line), indicating that the superposition 
A/± (|a) ± I — oi)) is completely interchanged between the systems, differently from the coherent state which does not 
swap to oscillator 1. 

Therefore, from Fig. 7(a) we conclude that the superposition state M± (|a) ± | — a)) prepared in a bad-quality 
system and protected from decoherence by coupling this system to a good-quality one, can easily be recovered 
in system 1 by switching off the coupling at the recurrence time = titt/A, n = 0,1,2,.... After this time, the 
superposition state will be in oscillator 1 with a fidelity less than unity due to the dissipative process. 

In Figs. 8(a and b) we plot the linear entropy for the joint state in Eq. (55) (<Si2 = 1— Tr/912) and the reduced 
states of oscillators 1 and 2 (<S^ = 1— Trp|). In these figures we employ the same parameters considered in Fig. 7, 
except that in Fig. 8(b) we use the ratio A/F = 2, instead of A/F = 20, to show clearly the dissipative dynamics. In 
Fig. 8(a) we analyze the recurrence-swap dynamics until around the correlation time tc, defined as the time when 
the entropy S\{S2) goes to about 0.1 at the swap (recurrence) time. In fact, as shown in Fig. 8(a), the minima of 
the entropy move away from zero due to the development of an inevitable correlation between the oscillators (due 
to the cross-decay channel) which thus become permanently entangled. This correlation time, estimated as the time 
when the minima of approach 0.1, is given by 

TC « ,| 1^ } + rr- (95) 

5|a|L£ Kit -It ) 

^From Eq. (95) we conclude that for the weak coupling regime, where 7/ = ij 1 the correlation time goes to 
infinity, i.e., the entropy Si{S2) always returns to zero in the swap (recurrence) time. Therefore, in the weak coupling 
regime the oscillators do not get permanently entangled (due to the absence of the cross-decay channel) and, since 
Si{S2) always returns to zero in the swap (recurrence) time, one is always able to recover a superposition state of 
a bad-quality oscillator coupled to a good-quality one. However, even in the strong coupling regime the correlation 
developed between the oscillators could not affect the process of recovering a superposition state of a bad-quality 
oscillator coupled to a good-quality one. Since the ratio tq / td obeys 

^ _ Hi +77) .Qf,X 

rn^ 5 E,(7/-77)' ^ ^ 

for Tel To ^ 1 one can always recover the superposition state M± (|a) ± | — a)), with a considerable fidelity, in spite 
of the process of entanglement between the oscillators. 

The overall pic;ture c;oniing from Fig. 8(a) is that the recurrence dynamics due to the strong coupling between the 
oscillators tends to restore the coherence of the initial states while the dissipative dynamics promotes the decoherence 
process. In Fig. 8(b) the linear entropy of the joint state 1S12, represented by the thick solid line, starts from zero, goes 
to a maximum due to the decoherence process and then returns to zero, since in the asymptotic limit both oscillators 
reach a pure state: the vacuum or some coherent state whose excitation depends on the amplification parameter F. 
Meanwhile, as shown in Fig. 8(a), the linear entropy St of the reduced state of oscillator I, oscillates between zero and 
0.5. The linear entropy 5i (^2 ) , solid (dashed) line, becomes zero when oscillator 1(2) assumes the state | —77)1(177)2), 
as can be computed from Eqs. (51a) and (51b). At the same time, S2{S\) bump into the thick solid line representing 
the linear entropy for the joint state of the system <Si2, from above, indicating that the superposition \t\)\) has swapped 
(recurred) to oscillator 2(1) on its way to decoherence. The maximal correlation between fields occurs at the points 
where the recurrence and swap curves cross, as illustrated by the dotted line representing the excess entropy, defined 
as 

J = 5i+52-5i2. (97) 

We also observe from the dotted line that the minima of the excess entropy X move away from zero, due to the 
development of an inevitable correlation between the oscillators. Returning to Fig. 8(b), after reaching its maximum 
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the correlation S12 decays exponentially as a result of the dissipation and the driving field, attaining zero in the 
asymptotic limit when there is no correlation between the fields described by stationary factorized coherent states. 
During the decay of correlation S12, the linear entropy Se of the reduced state of oscillator i does not attain the value 
»Si2 , due to the correlation described by the excess entropy I. 

IX. CONCLUSION 

We have presented a comprehensive treatment of the coherence dynamics in a network composed of two coupled 
dissipative oscillators. First, we have derived a master equation for both regimes of weak and strong coupling between 
the oscillators. In the weak coupling regime the dissipative mechanism of the individual oscillators is not significantly 
affected by their interaction, which appears only in the von Neumann term of the master equation. However, in 
the strong coupling regime the time evolution of the density operator of the joint system is modified by a cross- 
decay channel represented by a Liouville operator £12^12 accounting for the coupling between the oscillators. The 
appearance of this cross-decay channel leads to interesting properties of the coherence dynamics of strongly coupled 
oscillators. 

After the mathematical development, we first analyzed the state-swap and the state-recurrence dynamics, i.e., the 
probability that each oscillator returns to its initial state and the probability of state-swapping between the oscillators, 
respectively. In particular, we have analyzed these processes for the case where both oscillators, prepared in the joint 
state J\f± {\a) ± | — a))^ (g) |ry)2, present the same quality factor T. In the weak couphng regime the recurrence and 
swap processes are fully accomplished apart from the relaxation of the field states due to dissipation. On the way 
to the strong coupling regime the recurrence process remains unchanged while the swap dynamics is gradually lost 
as the coupling gets stronger. In the strong coupling regime the swap dynamics is lost due to phase mismatching 
between the coupling parameter A and the field-shifted frequencies we, as explained in Sec. V. Evidently, in the strong 
coupling regime the field states recur much more often than in the weak coupling regime before the relaxation takes 
place. 

Next, aware of the decoherence dynamics in the weak coupling regime, governed by the usual master equation (21) 
where a Liouville operator Cepi2 accounts for the effect of the reservoir on oscillator £, we turn to the decoherence 
process in the strong coupling regime. In this regime, the normal-mode frequencies {ujf^ ~ ujf zL X) of the coupled 
systems are strongly shifted, away from the oscillator frequencies wio, to regions of the reservoir frequency spaces 
where the spectral densities may be significantly different from that around iveo- As the spectral densities of the 
reservoirs play an important role in this regime, three different spectral functions were considered for our analysis of 
the decoherence process. When the normal-mode frequency uj'^ is shifted to regions around the origin of the frequency 
space, even for Markovian white noise the coupling parameter of oscillator £ with its reservoir, around lu^ , becomes 
half its value aroimd uj'^ (i.e., 7^ = r/4). When considering, instead, a Lorentzian spectral hmction, which approaches 
zero around the origin of the reservoir frequency space, the coupling constant 7^ becomes even smaller, resulting in 
decoherence times for the joint and the reduced system states significantly longer than the result computed for a 
single system plus reservoir. In Sec. VI the decoherence times for three different states of the composed system are 
computed, considering the three distinct reservoir spectral functions. Apart from the dependence of the decoherence 
time upon the reservoir spectral density, we stress that the coherence decay of the eigenstate of the normal mode u)^ 
becomes slower than that for the eigenstate of the normal mode w^. 

As discussed in Sec. Ill, apart from the possibility of considering particular physical systems with appropriate 
reservoir spectral densities, it is possible that specific spectral functions could be achieved through engineered reser- 
voirs. The difficult task of engineering strong interactions between the oscillators, together with the achievement of 
specific reservoir spectral functions, are the most sensitive problems in the way of the physical implementation of the 
network here proposed. However, this proposal might provide a motivation for future theoretical and experimental 
investigations. In particular, the proposal could be applied to test the Markovian white noise approximation or to 
probe the reservoir spectral functions. 

When considering the two systems to have different damping rates Fi and F2 (with A ^ Fi,F2) we demonstrated 
that the coupling between the systems, independently of its strength, makes the good-quality system worse and 
the bad-quality system better. This result can be employed to improve the quality factor of a cavity, and thus to 
protect quantum superposition states generated, via atom-field interaction, in open bad-quality cavities coupled to 
closed good-quality ones. Evidently, this results holds for the weak coupling regime because we have assumed that 
the coupling strength A is significantly larger than the system damping rates even for that regime. In fact, when 
considering A ^ F^, the computed improvement of the quality factor of a system follows from the recurrence- swap 
dynamics which take place many times before the relaxation time, protecting in a good-quality system a field state 
originally prepared in a bad-quality one. 
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Finally, wc have developed a careful analysis of the entropy excess in our network. Supposing that the joint system 
is prepared in the factorized state 1^12) = A/± (|ct) ± | — ® |?7)2, we observed that in the weak-coupling regime the 
recurrence-swap dynamics take place uninterruptedly until the relaxation process rules it out, i.e., the entropy excess 
always returns to zero in the recurrence and swap times. Therefore, in the weak-coupling regime both oscillators 
always get disentangled in the recurrence and swap times. Differently, in the strong coupling regime a correlation 
is developed between the states J\f± {\a) ± \ — a)) and \rj), in a such way that they get permanently entangled after 
a time interval we have called the correlation time. In this way, one cannot recuperate - with a fidelity equal to 
unity - a superposition state Af± (la) ± | — a)) prepared in a bad-quality system and protected from decoherence 
through a strong coupling of such system with a good-quality one. However, wc computed the amplitude \a\ where 
the correlation time becomes longer than the decoherence time, allowing the recuperation of the superposition state 
A/± (I a) ± I — a)) with a good fidelity. 

One of the main restrictions argued against quantum computation is that even if each individual logic unit were 
only slightly affected by the decoherence process, the coupling of a large number of such logic cells, for an actual 
implementation such as the factorization of large numbers [14,39], would decrease the decoherence time in such a 
way that the whole computing process would be seriously compromised [16]. Here we have shown, at least for two 
strongly coupled sites, that the decoherence process is not dictated simply by the excitation of the state involved in 
the logic operation and the damping decay rate of the logic cells. In a strongly interacting quantum network the 
decoherence process is delayed by the cross-decay channel arising from the strong coupling between the logic cells. In 
fact it may occur that for many coupled sites, as for the two coupled cavities presented in this paper, the decoherence 
time increases. Therefore, the discussion about strongly interacting oscillators we have presented here is central to 
understand the decoherence process in a quantum network; specifically, our conclusion is that the decoherence time 
depends not only on the excitation of the entanglement involved in the logic operation but also on the coupling 
strength between the logic units. 
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figures caption 

Fig. 1 Sketch of the coupled dissipative oscillators, with oscillator 2 submitted to a classical driving field. 

Fig. 2 Damping function 7^ (x) assuming a Lorentzian coupling Vg between oscillator £ and its respective reservoir. 
In the weak coupling regime the function 7£ (x) is centered around loi (dotted line). As A increases, the damping 
function splits into two Lorentzian functions whose peak heights are half the original value F^ (solid line). On the 
way to the strong coupling regime the two peaks can be clearly distinguished as shown by the dashed line. 

Fig. 3 Spectral density of the reservoir ag (x) for (a) Markovian white noise, (b) a Lorentzian spectral density, and 
(c) a wide Lorentzian spectral density. The system-reservoir couplings around oj^ are represent by the shaded regions. 

Fig. 4 State-swap probability Ps{t) (dotted line) and recurrence probability PR{t) (solid line) as a function of the 
scaled time Xt, for the factorized state IV'12) = A/± (|a) ± | — a))^ (g) 177)2, taking a = 77 = 1 as real parameters. In 
(a)-(c) relaxation and driving field are disregarded and we assume (a) weak coupling, (b) an intermediate coupling, 
and (c) strong coupling regimes. In (d) and (e) we set F = but dissipation is included, considering (d) weak coupling 
and (e) strong coupling regimes. Finally, in (f) dissipation is disregarded and the driving field turned on. 

Fig. 5 Decoherence dynamics for the factorized state \1jj12) = A/± (la) ± | — a))-^ $5' |?])2: assuming both oscillators 
with the same damping factor. Coherence decays of oscillator 1 (dashed line), oscillator 2 (dotted line), joint system 
(solid line) and an isolated oscillator with damping rate Fi (dashed-dotted line) are indicated. 

Fig. 6 Decoherence dynamics for the factorized state IV'12) = ■Af± (|q:) ± | — a))-^ ® 177)2, assuming that the damping 
rate of oscillator 1 is larger than that of oscillator 2. The dashed (dotted) line refers to the decoherence dynamics of 
oscillator 1(2), while the solid (dashed-dotted) line represents the coherence decay of the joint system (an isolated 
system with the damping rate of oscillator 1), respectively. The thick solid line represents the coherence decay of a 
superposition state A/± (jo;) ± | — a))2 prepared in an isolated oscillator 2. 

Fig. 7 Probability 'Pe{t) ~ Ipeit)] ''Pm), of finding the superposition state = J\f± {\a) ± | — q:))^ and the 
coherent state \ip2) = \r])2 in oscillator £. We consider (a) £ = 1 and (b) £ = 2. In these figures the strong coupling 
regime is assumed together with identical dissipative systems and Markovian white noise. 

Fig. 8 Linear entropy for the joint state (1S12 = 1— Tr/512); represented by the thick solid line, and the rodTiccd states 
of oscillators 1 and 2 {Sg = 1— Trp^), represented by solid and dashed lines, respectively. The dotted line indicates 
the excess entropy. In Fig. 8(a) the recurrence- swap dynamics is plotted until around the correlation time tc when a 
permanent correlation is developed between the oscillators. In Fig. 8(b) we are concerned with the relaxation of the 
network. 
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